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Abstract 

We study the and v T decays in the early epoch of the universe. 

If m VT > 2m e , there would be a CP asymmetry between v T — > e + + 

e~ + u e and u T — > e + + e~ + v e . The resultant CP non-invariance 

is a function of temperature and density and can reach 10 for a 

reasonable temperature range, but it is noticed that if m UT > 2m^ ~ 

200 MeV, the CP non-invariance can be much larger. 
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I. Introduction. 

As a possible candidate of the hot dark matter, the characteristics of neutrinos 
need to obey many constraints from both cosmology and experiments on the Earth. 
The recent research indicates that the upper limit of the dark matter neutrinos must 
be less than 10 eV Q. The earth experiments give an upper bound on v e mass as 
m Ve < 7.2 eV |§, the measurements by a PSI group |§ has given an upper limit of 

as m v < 160 keV. For z/ T , there are several groups which have obtained various 
upper bounds as [||] in Table I. 
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Method 


1-D 


1-D 


2-D 


2-D 


m„ e (<95% CL) 


32.6 


31 


74 


23.8 



Table I. The mass is in unit of MeV. 



The stringent constraints demand that either the neutrino as the hot dark matter 
component is lighter than 10 eV, or can decay fast enough, so that its relic contents 
are not substantial in the present universe. Of course, there is possibility for ex- 
istence of heavy unstable neutrinos to which the cosmology sets lower bounds ||. 
But a crucial problem is to investigate the decays of the lighter neutrinos and v T . 

The cosmology demands that corresponding to a few MeV, the minimum lifetime 
is r < 2.5 x 10 8 sec. for Dirac neutrinos and r < 4.3 x 10 8 sec. for Majorana neutrinos 

As long as v T is heavier than 2m e ~ 1.0 MeV, it can decay via v T — > v e + e + + e~ 
and the decay rate at tree level can be written as 
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where the electroweak (EW) corrections and the mass effects of m e are omitted, U Ut£ , 
U Uee are the Cabibbo-Kobayashi-Maskawa matrix elements. Provided vtl Vt ~ 31 MeV 
and U VTe ~ 0.01, U Vee ~ 1, the lifetime r can be about 11 sec. Therefore if m„ T is 
sufficiently heavy, it can satisfy the cosmology constraint. However, if m Vr were 30 
times lighter (say, 1~2 MeV), its lifetime would be close to the cosmology bound. 

A more serious problem exists for if its mass is about 160 KeV. Due to the 
phase space restriction, cannot decay via any reactions at tree level and the most 
probable decay mode is z/ M — > z/ e 7 which can only occur via the penguin diagram. 
Because of the loop suppression and the GIM mechanism, the decay rate would be 
very small, if the temperature effects are not taken into account. In fact, the lifetime 
can be as long as 10 22 sec. i.e. 3 x 10 14 yrs. 

Nieves et al. systematically studied the properties of neutrinos in medium with 
finite temperature and density J7|, especially, they investigated v T — > z/ e +7 and found 
a tremendous change, because then the GIM suppression is dismissed, so that the 
decay rate can be nine orders larger than that at zero temperature [§] . Following the 
trend, we will employ their results to analyze the lifetime, provided ~ 160 KeV. 
We further study the temperature and density effects on the decay v T — > v e + e + + e~ . 
Obviously, the temperature and density effects occur at the propagators which only 
exist in the penguin loop. At zero temperature, the loop correction is small with a 
factor a /it ~ 2 x 10~ 3 with respect to the tree level. For higher temperature and 
density which exist at the early epoch of the universe, there are two more energy 
scales besides the mass of z/ T , e.g. the temperature T and the electron chemical 
potential fi which corresponds to the density. Therefore once T, fi > m^, there is 
a chance to elevate the amplitude by orders and this enhancement may compensate 
the loop suppression. So we investigate the temperature and density effects via 
loops, concretely the penguin diagrams. 

The penguin diagram at zero-temperature has been carefully studied by many 



authors 0. Nieves et al. studied the penguin contributions to v T — > v e + 7, where 
7 is a real photon satisfying on-shell condition q 2 = 0. On contrast, in the process 
v T — > u e + e + + e~, the propagating photon is a virtual intermediate boson and 
q 2 7^ 0, therefore the results would be expected to deviate from that in v T — > 
v e + 7. Moreover, it is well known that as q 2 7^ 0, the penguin loop may possess 
an imaginary part which induces a CP asymmetry for v T — > v e + e + + e~ and 
v T — > z/ e + e + + e~ 0]. The temperature and density effects could modify the 
CP non-invariant phenomenon. In fact, any CP non-invariant source in the early 
universe is interesting and worth careful studies. 

Our paper is organized as following. After the introduction, we give the neces- 
sary formulations, in Sec. Ill, we present our numerical results, while the last section 
is devoted to a brief discussion and conclusion. 



II Formulations 



(i) For u e + 7. 

Because the final phase space forbids the reaction — > v e + e + + e~ as long 
as < 1 MeV, so the most reasonable decay mode for is the radiative decay 
*4 + 7. 

Nieves et al. [0 || assumed that only electrons (positrons) in the atmosphere, 
for the case of a nonrelativistic (NR) electron background, 



T (NR) 



1.39 x iy rF(V) 



1Q 24 cm -3 



(-)' 

V m / 



(2) 



whereas for an extreme relativistic (ER) background, 

p(ER) , rp s 4 

= 1.5 x 10 9 rF(V) — 

T \m, 



(3) 



where r = \U* Ve U eu \/\U* Ue U TV \ which is a large number about 10 4 ~ 10 5 unlike that 
in v T — > z/ e 7 it is close to unity, T is the width at zero temperature. 



F(V) = (1 -V 2 ) 1/2 



V VI- V 



where V is the three- velocity of the decaying neutrino and F(V) can take values 
between 1 and 1.55. The superscripts in eqs.(0, |3|) denote the non-relativistic and 
extra-relativistic cases respectively and n e is the electron density in the surrounding, 



m is the mass of u^. 



The zero temperature T is expressed as 
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Obviously due to the GIM mechanism, only r contributes substantially at zero 
temperature. 

(ii) For v T — > v e + e + + e~. 

The dispersive part of the amplitude comes from both the tree level and the 
the penguin diagrams, while the absorptive part only comes from the penguin. It 
is noted that besides the absorptive phase determined by a ratio of the absorptive 
part and the dispersive part, there are the CKM phases and they would be dealt 
with separately. 

The propagator of fermions at finite temperature and density can be written as 

[11] 

1 



Sf(p) = (j> + m) 
where 



p2 _ m 2 _|_ i t 



+ 2irS(p — m )r](p ■ v) 



(5) 



9{x) 6{-x) 

and (3 = 1/kT, \x is the chemical potential and v a is the four- velocity of the medium, 
in the Lab frame v a = (1,0). 

The Feynman diagrams are shown in Fig.l where (a) corresponds to the tree 
diagram, obviously at this energy scale, (c) gives rise of negligible contributions only. 
In the real photon emission case where q 2 = 0, 5(p 2 — m 2 e ) and 5({p — g) 2 — m 2 ) cannot 
be non-zero simultaneously as long as m 2 ^ 0, therefore, in the reaction v T — > v e + 7 
the penguin diagram does not possess an absorptive part. However, as q 2 7^ 0, 



product of the two S— functions can be non-zero, therefore an absorptive part of the 
penguin diagram emerges. For the energy scale of m Vr ~ 31 MeV, only electron-loop 
contributes an absorptive part, but not /i and r, because m UT < 2m^. It is a well- 
known fact in the zero-temperature field theory, but for the finite temperature and 
density situation, their effects on the absorptive part appear via rj(p ■ v)t]((p — q) ■ v). 
The calculations for the absorptive part concern an integration 

d* P 5(p 2 -m 2 M(p-q) 2 -m 2 e )f(p). (7) 



which was investigated by Kobes and Semenoff ]TT] . 

Thus, the total temperature and density dependent amplitude would be 

M = -i (^=j w^.7«(l - 75K T e7 a (l - 75 )e — _^ M ^ U e^ U eu r 

+ (-z)^ e r^Iey*e, (8) 
T 

where r^ 6 -^ stands for the loop contributions of e, /x and r as 

r< B//) = E K eff) + ^ eff) )Ku^ 0) 

l=e,n,T 

and the superscripts D and A represent the dispersive and absorptive parts of the 
penguin diagram respectively. 

In the derivations, we follow Nieves et al. [7] to use their conventions and de- 
composition as 

T a(3 = M T R af3 + M L Q a/3 + M P P af3 , (10) 



where 



Ra/3 = ga/3 ~ Qa/3, (H) 

ft,^ (12) 

Pa/3 = -Q^^a/3\fiu ( l^V U , (13) 



and g a p = g a p - v a = g a pv 
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Then the following calculations are straightforward and standard, even though 
tedious, we ignore the details in the work to save space. 

Thus the amplitude squared can be easily written down in a standard framework. 
The decay width is 



-itotal 



d 3 p Vr 1 r d 3 p e + 1 d 3 p e - 1 d 3 p Ve 1 



(2ir) 3 Ty T ' 4E„ T J (2n) 3 2E e+ (2tt) 3 2E e - {2n) 3 2E Ue 
■\M\ 2 5\p UT -p e+ - Pe - -p Ue )(l - /„J(1 - / e+ )(l - / e -), (14) 

where /j's are the distribution functions of fermions, (1 — fi) denotes the Pauli 
blocking of the produced e + , e~ and v e . It is suggested that in the early universe, 
one may use y/—g where g is the determinant of the metric matrix, instead of , 



2EV y ""^^ y ^ ^^^^u ^ "^«v,^ ^ 2£ . 

but here we assume that at this concerned time period of the universe evolution, 
the universe is flat enough and •>/— p is close to unity, 
(hi) The CP asymmetry. 

For v T — > e + + e~ + z/ e , the amplitude can be recast into another form which is 
more appropriate for the CP asymmetry discussions as 

M = A{ 0) e i5 * + A[ e) e i5 *e i9 ° + A M e tS » + A[ r) e i5 % (15) 

where 



A{ 0) = \A (0) u: u U eUT \, (16) 

(17) 



AD)2 A (A)2 jr JT 



A '(V) = \ A M U Ue U ^r\, (18) 

A[ T) = \A {T) U* TV U TVT \, (19) 

and A( ) is the tree level quantity, while Aq^s correspond to the loop as I — e, /i,t 
respectively. 

The absorptive phase 9 e which is equivalent to the strong phase in K — > 2n 
process, comes from the absorptive part of the penguin loop as 

tarr 1 ^-^, 9 li = 9 T = 0. (20) 

Ae 



Since m Vr < 2m^, only the electron loop contributes an absorptive part according 
to the Cutkosky cutting rule [[12] and the weak phases S e , 5 M , 5 T originate from the 
CKM phase. Under a CP transformation, the weak phases change signs, while 9 e 
does not. The interferences give rise to a CP asymmetry. 
The CP asymmetry is proportional to 

JUk(dLIPS k )(\M\ 2 -\W) 
a Jllk(dLIPS k )(\M\* + \M\*) 

_ J UkjdLIPSk)^ 2|^| |4| sin(^ - Sj) smjOi - Bj) 

" /n fc ^/^)[E i j[|^| 2 + 2|^||^|cos(5 J -^)cos(ft t -^)]' 1 J 

where dLIPSfc denotes the Lorentz Invariant Phase Space of final particle k. 

The numerator comes from interferences among all the terms (tree level and 
e, fi, t loops), obviously, it is zero unless both the weak and strong phases of the 
two terms are different. Therefore in our case, the numerator is proportional to 

2|A / (e) |P' (M) | sin(5 e - <y + \A[ T) \ sm(8 e - 5 T )\ sm6 e . (22) 

It is observed that there is no contribution to the CP asymmetry from the in- 
terference between the tree amplitude and A\ because they either have the same 
weak phase (with A', e \) or the same absorptive phase (with A',-, and A',^,). As some 
authors suggested that in the surrounding there were no \i and r leptons, then ^4',^ 
and M r \ are the zero temperature amplitudes which are somewhat smaller than that 
with higher temperature and density. 

We can argue that since m T ~ 1.7 GeV ^> m M , one can expect it to decay very 
fast, so that there is no r— lepton in the atmosphere, but probably a certain amount 
of muons. In our later calculations we assume that both \x and r exist when the 
reaction takes place. 

III. The numerical results. 

In the calculations, we study the high temperature and density effects for the 
early universe. The studies on the evolution of the universe suggested that the 
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reasonable temperature and density limits may be 4 x 10 K and /i = 600 MeV, so 



we take them as the extreme bounds \T3\. In the calculations we adopt the CKM 
matrix elements as given in the data book [13] and ref. JTi5fl . 

(i) For the decay, the zero temperature calculations predict a rather long 
lifetime as r v ~ 3.5 x 10 12 yrs which is longer than the lifetime of our universe. 

The calculations show that at T = 10 12 K, the lifetime of decreases to 2 x 10 3 
sec. But since the temperature of the universe cannot hold so high for more than 
a few seconds, thus one can expect that an average lifetime of would be 10 s 
sec, which roughly coincides with the cosmology requirements. We will discuss this 
consequence in the last section. 

(ii) For v T decays. 

There are two important decay modes v T — > e + + e~ + v e and v T — > v e + 7. At 
zero temperature, because the former reaction can be realized via a tree diagram, 
so its contribution dominates over the radiative decay mode which can only occur 
via loops. 

The finite temperature and density effects change the characteristics of the decay 
modes. The radiative decay was studied by Nieves et al. in detail, while we investi- 
gate the temperature and density effects on v T — > e + e~v e in this work. Because at 
high temperature and density, the new energy scale T and the chemical potential 
would replace or partially replace the mass term, so can result in a larger decay 
rate. Our numerical results indicate that as T and \i get large enough, the loop 
contribution can be one times larger than the tree contributions, (see Fig. 2). 

In fact we have also found the dependence of the decay rate on the density 
as T increases fast with the density increment. It is the main aim of this work. 
The dependence on the chemical potentials are shown in Fig. 2, the three curves 
correspond to various temperatures. Its significance would be discussed in next 
section. 

It is noted that the penguin contributions increase fast with increments of tem- 
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perature either, for // = 0, at about T = 6 x 10 11 K, the tree amplitude and the 
penguin have the same magnitude, then the penguin contribution takes over and 
dominates. 

(hi) The CP asymmetry. 

One could expect an absorptive part emerging at zero temperature as an absorp- 
tive angle is resulted |17| 



tan d P 



G F 2a 
V2 3 



(23) 



which is about -0.4, but as the temperature is not zero, this ratio becomes very 
small that 9 P ~ 10 -6 . Concrete values are shown in Table II. 
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a(xl0~ 7 


4.1 


4.1 


4.08 


3.98 


3.2 


2.6 


2.1 


1.5 



Table II. in the table, T = 4 x 10 12 K, \i is in MeV. 



In fact, in the interaction 

/ dW)5((p - q f) eP{E l ) + ieP{E l ) + i {2,fF{ P \ (24) 

where we ignore the small electron mass and F(p) is a function of p, the distributions 
are generally smaller than unity, on contrast, at the zero temperature, it is 

j d% P 2 )5{{ P -qf){nfF{ P ), (25) 

so it can be expected that the finite temperature and density effects cannot make 
the integration values to be much larger than the zero-temperature case. Especially, 
as the chemical potential /i gets larger, the second distribution in ( |24|) would de- 
crease very steeply, thus the integration turns smaller accordingly and so does tan# e . 
Therefore, we should think that the absorptive part remains at the same order as 
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the zero temperature value. Our careful numerical computation confirms it. Thus 
a rapid increase of the dispersive part of the penguin loop with temperature and 
density causes tan# e to decrease fast and the effects disfavor the CP asymmetry. 
We will discuss it in next section. 

IV. Discussion and conclusion. 

We have evaluated the temperature and density effects on the and v T lifetimes 
and estimated a possible CP asymmetry a for various T and \i. 

(i) For the lifetime, it can only decay via the radiative mode, we estimate its 
lifetime as about 10 3 sec. as m v = 160 KeV for very high temperature. However, 
obviously, such high temperature cannot hold for so long and would decrease very 
fast, and so does the decay rate since it is proportional to T 4 . Considering this, 
lifetime is close to 10 8 sec. which is the bound set by cosmology. It is possible that 
as < 160 KeV, its lifetime bound would be even larger, then one may conclude 
that to be consistent with the cosmology, either the m v is much lower than the 
upper bound to be close to 10 eV, or there is some new physics which speeds the 
decay, for example axion exists etc. 

(ii) The temperature and density effects on the v T lifetime. The process v T — > u e + 
7 was carefully studied by Nieves et al. ||. We study the effects on u T — > e + e~v e 
for which there is a tree diagram. For v T — > e + e~u e , the temperature and density 
effects occur via the penguin loop which is 10~ 3 — 10~ 4 orders smaller than the tree 
amplitude at zero temperature. The numerical results indicate that the dispersive 
loop contribution is enhanced to the same order as the tree amplitude. Moreover, 
at T = 10 12 K, the two decay modes are of the same orders 

We have found that as temperature increases, the rates of both modes increase, 
below a certain temperature (about T = 10 12 K) the rate of v T — > e + e~v e is larger 
than that of v T — > v e + 7, at the temperature, they are at the same order, but 
while the temperature further increases the radiative decay has a larger rate. The 
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concrete temperature also depends on the density, but the trend is obvious. As Babu 



et al. pointed out that the window for MeV neutrino might be closed |TH] , if so as 



aforementioned, m„ T gets to keV order, its lifetime would become as troublesome as 
for z/ M of 160 keV. 

(hi) The CP asymmetry. 

At zero temperature, it is known that the penguin diagram induces a CP violation 
which can be observed at B- and D-decays. As estimated, the absorptive angle which 
is equivalent to the strong scattering phase, is relatively large. It is natural to ask 
if it can result in a CP non-invariance at finite temperature and density which exist 
in the early universe. 

Since v T — > z/ e + 7 does not possess an absorptive phase because q 2 = 0, it cannot 
induce a CP violation. Whereas, for v T — > e + e~u e , there indeed exists an absorptive 
angle, but there is the tree amplitude which is overwhelming at zero temperature. 
Therefore, even though the absorptive angle is large, the CP violation is very small 
because the tree contribution exists at the denominator in expression (pl|). When 
the temperature is very high, the penguin contributions are much enhanced and 
can reach the same order as the tree contribution, however, it is observed that only 
the dispersive part increases, while the absorptive part remains at the same order. 
Thus the absorptive angle 9 e turns out to be rather small and the total CP viola- 
tion at high temperature and density keeps the same order as at zero temperature 
and remains to be small. Therefore, even though the decay can serve as one of the 
CP non-invariance sources, but definitely is not the main one to explain the matter 
dominance over anti-matter. 

Acknowledgment We would like to thank prof. R.K. Su for helpful discussions. 



Figure Captions 
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Fig.l The Feynman diagrams for v T — > e + + e + z/ e , where (a) is the tree diagram. 
Fig. 2 The dependence of the decay rate on the chemical potential, the three curves 
correspond to T = 1 x 10 10 K (solid), T = 5 x 10 n K (dashed), and T = 4 x 10 12 K 
(dotted), respectively. 
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